Analytic Relations between Localizable Entanglement and String Correlations 

in Spin Systems 
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We study the relation between the recently defined localizable entanglement and generalized cor- 
relations in quantum spin systems. Differently from the current belief, the localizable entanglement 
is always given by the average of a generalized string. Using symmetry arguments we show that 
in most spin- 1/2 and spin-1 systems the localizable entanglement reduces to the spin-spin or string 
correlations, respectively. We prove that a general class of spin-1 systems, which includes the Heisen- 
berg model, can be used as perfect quantum channel. These conclusions are obtained in analytic 
form and confirm some results found previously on numerical grounds. 
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Entanglement is one of the most specific features 
of Quantum Mechanics (QM), first recognized by 
Schrodinger , that has been realized in recent times to 
open the way to some of the most important applications 
of QM like teleportation, dense coding, quantum cryptog- 
raphy and computation Since quantum information 
is based on qbits and, more generally on finite-level quan- 
tum systems, spin chains emerge as natural candidates 
for quantum devices. Moreover, it has also been sug- 
gested recently that the quantum information perspec- 
tive could give a better understanding of the structure 
of the ground state and of possible quantum phase tran- 
sitions in Condensed Matter systems 0, 0, 01 hi this 
context Verstraete et al. have introduced the notion of 
localizable entanglement (LE) [f|, which is a very at- 
tracting new concept for its possible applications to the 
study of quantum phase transition. In addition, the LE 
has a precise physical meaning to the extent that it can 
measure the performance of certain kinds of quantum re- 
peaters 0,0. 

The LE of an interacting multipartite system is the 
maximum amount of entanglement that can be localized 
on two particles by performing local measurements on 
the remaining (assisting) particles. Here we stick to the 
class of projective Von Neumann measurements. The LE 
on the pair i, j of an N + 2 particle state \ip) is defined 
as 



while the upper bound is the concurrence of assistance 



L i,j = nmxyV.E (\<p a )) 
{\ s )\ „ 



(1) 



The maximum is taken over all local measurement bases 
{|s)} of the N assisting particles, \ip s ) is the normalized 
two-qbit state at sites i and j obtained after performing 
the measurement \s) with outcome probability p s . As a 
measure of entanglement E (\<p s )) for two-qbit states \tp s ) 
we use here the concurrence, first defined in [9J- Find- 
ing the optimal basis which satisfies Eq. (JTJ is in gen- 
eral a formidable task. However in upper and lower 
bounds have been derived, the lower bound being given 
by the maximal connected correlation between sites i and 



Interestingly, it has been observed numerically that in 
most spin- 1/2 systems the LE equals the lower bound 
. Similar numerical studies have shown [ill ll2T ] that 
the LE assumes its maximal value 1 for the ground state 
of the spin-1 Heisenberg model with spin 1/2 at the end- 
points. The same result has been proved exactly Il2l 
for the AfHeck-Kennedy-Lieb-Tasaki (AKLT) model [laf . 
Instead the 0-deformed AKLT model ^3] exhibits an 
exponentially decaying LE with a finite entanglement 
length £e, which is defined in analogy with the correla- 
tion length. For the latter model the LE has been consid- 
ered an essential indicator of the transition taking place 
at (j) = 0, since here £e diverges while the standard corre- 
lation length is smooth. To rule out the possibility that 
the LE can be ascribed to the string order correlation, 
the authors of Ref. have given an example of a state 
for which £g diverges even in absence of string order. 

In this Letter, it is shown that the LE is given by a 
generalized string whose form depends on the symmetry 
properties of the state. In particular the following ana- 
lytical results are given: i) in spin 1/2 systems we give 
arguments that show the tightness of the lower bound 
ii) for a large class of spin 1 systems, that include the 
Heisenberg model, the LE is proved to be 1, i.e. maxi- 
mal; iii) in the (^-deformed AKLT model we show that 
the LE is connected with a string order correlation. 

Throughout the text we consider a system of N as- 
sisting particles labeled from 1 to AT and two qbits on 
which we localize the entanglement that, without loss of 
generality, we place at the endpoints 0, N + 1. 

Once an optimal basis {|s) = \sx . . . sjv)} is found, the 
LE of an + 2 particle state \ip) is given by 

L(|V>>) = El P (l<M)l> \<f>.) = (s\1>), 

s 

where the preconcurrence P of an unnormalized two-qbit 
state \<j>), is defined by P (\<f>)) = (<p*\a y ® a v \(f)), and 
the complex conjugate is taken in the standard basis. It 
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is then clear that the LE can always be written as an 
expectation value 

Lm = M°%GMo» N+1 \r) (2) 

where the operator is given by 

G( s )=^|.s)( S *|sign(P(|0 s ))), (3) 

s 

and the sign function is sign (z) = zj \z\ for nonzero z, 
sign (0) = 0. As it is evident from the definition, the G^ 
operator depends on the optimal basis as well as on 
However, in most cases G^ s ' takes a manageable form 
that depends on the symmetry property of the state. On 
top of that, if the signs in Eq. J2J factorize into local 
terms, then G^ 8 ' becomes a string GW — G\Gi ■ ■ ■ Gn- 

In order to find the optimal basis {|s)} we take an in- 
finitesimal variation over all possible local unitary trans- 
formations obtaining the following set of extremal equa- 
tions: 

^Im[rf( S )sign(P(|^)))] = (4) 

8 

rf(s)=0, forP(|0 s » = O, (5) 

where we have defined 

Tf (s) = «K K\s){s*\ + \s){s*\ (rf )*] < +1 |V*>, 

and the r M 's are the generators of the unitary single- 
site group SU (D) (D is the dimension of the single site 
Hilbert space), /i = 1, . . . , D 2 — 1 and i = 1, . . . , N. 

In what follows we shall discuss separately spin 1/2, 
spin 1 and matrix product states. 

Spin 1/2. In this case the r^'s are the Pauli matrices. 
On the basis of Eqs. (@J and (0, it is straightforward to 
show that if is real and invariant under it rotations 
about the a-axis, P a = P3 i= Jj ) 1 cr" , a — x,y,z then an 
extremal is provided by the a a basis. This result can 
readily be applied, for example, to the Ising model in 
transverse field and to the XXZ Heisenberg model. 

Let us consider first the Ising model in transverse field: 

H = -AE«+i-E^- ( 6 ) 

i i 

For finite N the ground state is unique for any A and 
belongs to the P z = 1 sector. For |A| > 1 the gap between 
the lowest state in the sector P z = — 1 and the ground 
state vanishes exponentially with N and hence becomes 
degenerate with it in the thermodynamic limit. As the 
basis of a z is extremal, we calculate now G {z \ For A > 
a Marshall's sign theorem ensures that, in the o z basis, all 
the coefficients of in the same parity sector are non- 
zero and have the same sign. Then the preconcurrence 
P (\4> s )) is positive (negative) according to whether the 



parity of \s) is —1 (1). The case A < is reduced to the 
previous one by transforming the Hamiltonian © with a 
7r rotation about the z axis on the even sites. Collecting 
results, for any non-zero A, we get 

N 

G^-sign(Af +1 ]>f. (7) 

This implies at once that the LE calculated in this basis 
is = sign (A) JV+1 (^|<Tg(T^ r+1 |V'), i.e. the max- 

imal classical correlation, that is the lower bound. For 
|A| > 1, in the thermodynamic limit we must be cau- 
tioned that the LE depends on the particular combina- 
tion of the twofold degenerate ground states one chooses. 
Let us discuss now the XXZ model 

H = E + °t^+i + A ^+i) ■ (8) 

i 

The Hamiltonian commutes with all P x ^^ z discrete rota- 
tions, meaning that a x, s , a y, s and a z, s are all extremal 
bases. Which one yields the highest LE must be assessed 
by explicit calculation. Exploiting the Marshall's sign 
property we get G^ = — (—1) N, P® N , where V\ N is the 
projector onto the subspace of zero magnetization for the 
assisting spins. This means that the LE assumes the form 

i (z) (l^» = -(-l) JV (^l« + il^). (9) 

To evaluate G^ we need the Marshall's sign property 
in the a x, s basis. This can be achieved with a local 
unitary transformation which depends on the value of 
A. For A > 1 the required rotation is ]Jj exp {i-Ku^ji) 
which acts only on one sublattice. Instead for —1 < 
A < 1 we must further apply a uniform rotation of 7r/2, 
Y\- exp (iTr(Jj 1 4) . For N is even, the result is 

N 

G (x > = -(-i)" /2 n°i> £ (x) = -<«<^+iM> a°) 

i=i 

for A > 1, and 

G^) = -1, L (x) = -(tl)\a%a v N+1 \il)), (11) 

for |A| < 1. Clearly = L^ x > for symmetry. 

In the critical regime, — 1 < A < 1, it is well known 
that the transverse correlations dominate at least asymp- 
totically, while the z—z correlations are the highest in the 
antiferromagnetic phase A > 1. The conclusion is that 
for the ground state of the XXZ model the a x 's basis al- 
ways yields the highest LE. As for the Ising model, the 
lower bound is attained, in agreement with the numerical 
simulations in Ref. |ll|. indicating that our local maxi- 
mum is indeed a global one. It is remarkable that on the 
each side of the isotropic point A — 1 the LE at finite size 
is given by different classical correlators. In particular it 
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has been observed in finite size systems that the LE be- 
tween nearest-neighbour sites displays a singularity [llj |. 
In our opinion, this fact is due to the different transforma- 
tions required by the Marshall's theorem and cannot be 
attributed to the Berezinskii-Kosterlitz-Thouless (BKT) 
transition occurring at the same point. 

Work is in progress to extend the method outlined here 
in presence of a magnetic field that breaks the P x sym- 
metry Nonetheless, the numerical computation of 
anticipates that for some values of the field, the LE is 
strictly greater than its lower bound. 

Spin 1. Here we consider spin one assisting particles 
Sf : i = 1, . . . , N, and two qbits at the endpoints. The 
extremal equations are still given by Eqs. 10} and (JSJ) 
where the r^'s are now the generators of SU (3), i.e. the 
eight Gell-Mann matrices . When both the state \tp) 
and the basis |s) are real, it is enough to consider in 
the extremal Eq. only the restricted set of purely 
imaginary generators of SU (3): (r p )„ „ = — ie M!yp , where 
t iivp - 1S ^ e Levi-Civita symbol. 

Given the richer structure of SU (3) with respect to 
SU (2), the invariance of the state through a rotation of 
7r around a-axis, 11" = e lwStot , is not sufficient to prove 
that the basis of S a is extremal, at variance with the 
spin 1/2 case. Instead, if the state \ip) is symmetric with 
respect to rotations of ir about two (and hence any) axes, 
using the following (anti)commutation rules 



0. 



Eq. (0} is satisfied provided that the measurement basis 
fulfills 



\s}(s\ 



AttS x 



= 0. 



(12) 



The solution of (|12J) is (apart from unimportant phase 
factors) {|5» = {|0>, \±)} = {|0>, (| + 1) ± | - 1)) /^}, 
since the matrices e mS ' are all diagonal in this basis. 

Surprisingly, without making any further assumption 
on the state it is possible to calculate exactly the 
LE in the basis {|s)}, finding that it reaches its maximal 
value, i.e. 1. The proof goes as follows. We write the 
state lib) as 



,a\siS2 ■ ■ ■ SN}\o-0O-N+l), 



(13) 



where \si) € {|0),|±)} and |er) is the a z, s basis for 
the two qbits. Due to the invariance of the state 
with respect to rotations of 7r about any axis, one has 
e i7r,St ° t V> = Vai'-, Pa — ±1 and the coefficients in Ijl3|l sat- 
isfy C s ,a = -PzOWjv+i (-l)" +l ^ +Tl -( s) C s ,o and C s ,a = 
—p x (—\) n ~^ +no ^C s ^ s . Here n 7 (s) is the number of 
times Si — 7 in the state \s1S2 ■ ■ • sn) and |cx)is obtained 



by spin flipping the state | a). The LE can be written as 
W = 2 £ \C s , n C Sill \+2j2 \Cs, U C s , u I = 

s s 



Now we use \C s , a C sra \ = - (_i)«+M+"°(s) C s ^C s - a and 



(-1) 



i + (s)+n (s) _ / I T-fiV iTrS; 



s\U"=i e 1 l s >- Finally, one obtains 



: Pa( 



x,y,z. 



(14) 



This implies that Eq. (0 is automatically satisfied. 
In this derivation it has not been necessary to com- 
pute the G operator which anyway is given by — 
—PyVYiiLi e llrS iV, where V projects out the states for 
which the preconcurrence is zero. Along the same rea- 
soning it can be shown that Eq. I|14|) remains valid for 
general complex 

The finding of Eq. (|14fl demonstrates that, in a spin 
1 state the entanglement is fully localizable on the end- 
points, with the only requirement of invariance with re- 
spect to 7r rotation around any axis. This property is 
shared by ground states of a large class of systems which 
can be regarded as perfect quantum channels. 

An example is given by the XXZ 5=1 Heisenberg 
model with single-ion anisotropy: 



H = 



SfS; 



i+l 



QV qV 

J; J; 



\S^ S r 



i+l 



d {sty 



(15) 

This model exhibits a very rich ground state phase dia- 
gram (see 0, ^| and references therein) with six differ- 
ent phases. According to numerical studies 0, the LE 
is one at the isotropic Heisenberg point. In the present 
work it is proved that the LE is always one in the entire 
region of parameters, in spite of the complexity of the 
phase diagram. Of course this makes unpractical the use 
of the LE for the detection of quantum phase transition 
in contrast to what suggested by some authors 0, 0] . 

One may wonder whether the maximal localizability 
of the entanglement is due to the presence of the two 
qbits at the borders. This occurs, for example, in the 
large-D region where the asymptotic ground state is 

(| T,0 0,4.) — I |,0,...,0,t))/V2for which the LE is 

one, while in absence of the spin 1/2 at the endpoints, the 
ground state tends to the unentangled state |0, . . . , 0). 

A little care is needed in the broken symmetry phases 
i.e. ferro- and antiferro-magnetic regions where the 
ground states are twofold degenerate. However the re- 
quirement of invariance under J\ x,v,z fixes the combi- 
nation of ground states in a similar manner as in the 
Greenberger-Horne-Zeilinger (GHZ) state with maximal 
LE. 

Matrix Product States. Matrix product states (MPS) 
are very interesting in many respects, in particular the 
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fixed point of the density matrix renormalization group 
algorithm (DMRG) yields an MPS state 0- As was 
shown in Ref. [12J, it is possible to calculate exactly the 
LE for any MPS state. This opens up the possibility 
of calculating approximately the LE for any one dimen- 
sional ground state. Here we revisit the examples intro- 
duced in in the framework of our formalism. This will 
allow us to clarify the link between LE and correlations. 

Using the convention of [12j we write a general MPS 
state with open boundary condition as 

\F A ) = C J2 \Pi — 0N)(l®At 3 «~.A*)\*-)o tN+1 , 

Pi,- ,Pn 

(16) 

where the A^'s are D 2 x 2 matrices which depend on 
the local basis C is a normalization constant and 
I*~)o,jv+i = (I U)o,jv+i - I 4t)o,;v+i) A/2 is the two- 
qbit singlet at the endpoints. 

The optimal basis |sj) can be found via the method 
shown in |l2| and, since in this case the preconcurrence 
takes the form - det (A SN ■ ■ ■ A S2 A Sl ), the G operator 
factorizes into local terms reducing to 



JjGi, Gi = J2^H^HA St ))\si){s*\ 



Hence the localizable entanglement is 

L{\F)) = -{F\alG x G 2 ---G N a y N+l \F 



(17) 



A well known example of MPS state is the ground 
state of the spin-1 AKLT model, characterized by the 
presence of a gap, exponentially decaying spin-spin cor- 
relation function and hidden topological (string) order 
[l3| . As any MPS, the AKLT ground state can be alter- 
natively written in the valence bond solid (VBS) formal- 
ism. In this representation, every spin one at position 
i is replaced by a pair of qbits i,i each one forming a 
singlet with its nearest neighbor, respectively i — 1 and 
i + The on-site couple of qbits is then projected back 
onto the spin-1 Hilbert space by means of a single 3x4 
matrix A, related to the 2x2 MPS matrices A 13 via 
(ft|A=(*-|^ft®l ? . 

The </>-deformed generalization \V^) 12] of the AKLT 
state is given by the matrix 



A = 



/V 





v/2 







ill 





The |V0=o) state is the rotation- invariant AKLT state, 
while a non zero (f> breaks the O (3) invariance down to 
O (2) invariance. As was shown in |l2j the optimal basis 
is {\s)} = {|0), |±)} for any <j>, which is the same as in the 
previously considered S = 1 case. It is straightforward 
to calculate the G operator and the LE which now reads 



N 
z=l 

(18) 

For cj> = the LE is equal to one, while for non zero <f> 
it decays exponentially with a finite entanglement length 
i E l = - lmiAr^oo In (L) /N. On the basis of Eq. O, the 
connection between the LE and the string order parame- 
ter (SOP) in the y direction becomes evident. Generally 
the (maximally extended) SOP's g a (N) are defined as 



N-l 



9 a (N) = -WIS? n 



S%\V<p 



i=1 



In the isotropic case = 0, all g x >y> z saturate to 4/9 re- 
vealing a breaking of a hidden Z2 x Z2 symmetry |l7j . 
On the other hand, for <p ^ the SOP in the z direc- 
tion tends exponentially to a non-zero value g z (00) = 

cosh (2(f>) + sj cosh (2(f) 2 + 3^j , while g» (and g x ) 

behaves as the LE decreasing exponentially to zero with 
correlation length £g. In the (j> = case the LE can be 
ascribed to the presence of string order. However in the 
deformed case, it was argued in that the LE is not 
connected to string order since the LE is short ranged 
while g z is non vanishing. Actually our arguments show 
that in this case a connection still exists between the LE 
and SOP, but in the y channel. 

An example of MPS state \rf) for which the LE is identi- 
cally one but all the SOP's are zero was conceived in 
and is defined as in (fTfi)l with A 1 = a z +a y , A 2 = a z — il, 
and is the standard basis which in this case coin- 

cides with the optimal one. It is straightforward to cal- 
culate the G operator for this state which turns out to 



be G = — (— 1) N 1 yielding a maximal LE: 



L(\v)) = -(-l) N (vKo« N+1 \f l *) = l. 



(19) 



It is clear that, in this case, the LE is not connected 
to string parameters but rather to a classical correlation 
that involves a complex conjugation. 

In conclusion, in this Letter we have considered the re- 
cently defined LE and its relationship with correlations in 
spin system. In quantum information theory the LE has 
a well defined physical meaning, while in condensed mat- 
ter it has attracted much attention for its possible uses 
to detect and characterize quantum phase transitions. In 
particular, we have shown that for spin-1 systems parity 
symmetry along the three axes is sufficient to assure that 
the LE reaches its maximal value one. This promotes 
spin-1 systems as realizations of perfect quantum chan- 
nels, but shows that the LE is insensitive to quantum 
phase transitions. Using similar symmetry arguments we 
have put in evidence that, in the spin-1/2 case, the LE 
equals the maximal correlation. Both of these findings 
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have been obtained analytically confirming the numerical 
results of Ref . • We believe that relaxing the symme- 
try requirements and/or analyzing higher spin systems 
opens the possibility of finding examples in which the 
LE does not reduce to known correlations. 
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